Nonlinear sigma model for optical media with linear absorption or gain 
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In the framework of the Keldysh technique, we formulate the nonhnear sigma model for disordered 
optical media with linear absorption or gain. The effective action for fluctuations of the matrix field 
about the saddle point acquires an extra term due to the nonconservative nature of the system. We 
determine the disorder-averaged Green-function correlator, which has a diffusion pole modified by a 
finite absorption/gain rate. The diffusion coefficient is found to be close to its value for conservative 
systems in the relevant range of parameters. In the medium with gain, the random-lasing threshold 
depends on the sample size. 
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I. INTRODUCTION 

The transport of waves through disordered matter has 
been a topic of recurring interest ever since the discov- 
ery of the Anderson localization in electronic systems [1] . 
Analogous phenomena have been subsequently studied 
for the transport of classical [2-5], matter [6, 7], and even 
seismic waves [8]. 

The research on classical-wave propagation in disor- 
dered media has been motivated by the conjectured pos- 
sibility of the localization of light. The results, such as 
the enhancement of dwell times due to resonant scatterers 
and, hence, lower energy-transport velocities [9] and the 
correction term in the Ward identity due to frequency- 
dependent scattering potentials [10], have shown that, 
while retaining many similarities, the behavior of light 
in disordered media differs from that of electrons in sev- 
eral important aspects. One of these aspects concerns the 
propagation of light in nonconservative disordered media. 
Such systems can be physically realized, for example, as 
random lasers [11-13], which have received much atten- 
tion recently. A promising research direction in this con- 
text are theories that combine description of wave prop- 
agation through disordered medium with the nonlinear 
laser equations [14-16]. 

The properties of light diffusion in absorbing media 
was studied using the photon transport equation [17- 
19]. In particular, it was argued that, in the parame- 
ter range of validity of the diffusion equation, the dif- 
fusion coefficient is close to its value in the conserva- 
tive medium. The treatment of light propagation start- 
ing from the wave equation has been mainly conducted 
via the self-consistent diagrammatic theory [20, 21]. In- 
teresting results, such as corrections to the bare diffu- 
sion coefficient due to the additional terms in the Ward 
identity [22, 23] and dynamics of Anderson localization 
in quasi-one-dimensional geometry [24] and open three- 
dimensional media [25] have been obtained by these 
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methods. An alternative description of classical wave 
propagation is provided by the so called effective models 
of disordered systems, commonly known as the nonlinear 
sigma model [26, 27] (NLSM). Being originally developed 
for electronic systems, the (supersymmetric) NLSM de- 
scribing light propagation in a conservative disordered 
medium was derived in Refs. [2, 28]. Later, the effects 
of an open boundary on the diffusion coefficient were 
studied [29] by using a similar model. Unlike the self- 
consistent theory of transport, effective models have not 
yet been applied to optical systems with absorption or 
gain. The effective models can be useful, e.g., in de- 
scribing special properties of light localization in such 
systems [30-32]. 



In the present work we formulate the Keldysh nonlin- 
ear sigma model [33] for the propagation of electromag- 
netic waves in nonconservative disordered media in the 
diffusive regime. Systems with absorption or gain are rel- 
atively simple to treat in the Keldysh formalism, which 
makes it possible to define an action needed for the field- 
theoretical description. By following the general scheme 
as outfined in Ref. [33] for electronic systems, we derive 
an effective NLSM action where we obtain a term due 
to nonconservativeness of the medium. A similar con- 
tribution was found in Ref. [29]; in that case the term 
originated from the openness of the system. 



Furthermore, by using the standard methods [33], we 
show that the light propagation can be described by a dif- 
fusion equation for nonconservative medium. The condi- 
tions under which the NLSM yields the diffusion equation 
are found to be equivalent to the restrictions imposed in 
the theory of transport equation [17]. Similarly, the dif- 
fusion coefficient that we derive is almost independent of 
the absorption or gain under these conditions. For the 
amplifying medium, we discuss the applicability of the 
linear-gain approximation and determine the threshold 
of random lasing. 
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II. KELDYSH APPROACH TO LIGHT 
PROPAGATION 

A. Partition function for nonconservative medium 

We consider an optical medium defined by a complex 

dielectric constant e(r, w). Restricting our theory to the 
TM modes in two dimensions, we describe the electric 
field by its normal component Ei^{r) — iujAi^ir) (in the 
Coulomb gauge), where ^^(r) is the normal component 
of the vector potential. We use the Gaussian units with 
the velocity of light c = 1. The transversality condition 
V • (eA) = with e varying in two dimensions, leads to 
the two-dimensional wave equation 



[V^ + e{r,uj)uj^]A^{r)=0. 



(1) 



For real e(r, w), this equation, and its complex conju- 
gate, can be obtained by setting to zero the functional 
derivatives 



SS 



— - n 
SA~ ' SA 







(2) 



of the action (Hamilton principal function) 



S[A,A*] 

= i/ dr^Hr,uj)oj'\A^{r)f - |VA.(r)|^] 
= ifc/ '^'■^^-^'■^ [e(r,u;)c.2 + V^] A^(r), (3) 

treating A and A* as independent functions. The ac- 
tion can be rewritten in the representation-free operator 
notation as 



(4) 



where the inverse Green function operator G^^ — 
e(r, co) uJ^ -\- in the (r, lo) representation and A (A'f) is 
the Hilbert-space vector ^^^(r) [A* (r)]. 

In order to construct the quantum Hamiltonian, one 
expresses the energy of the system in terms of the vector 
potential, ^^^(r) and ^j^(r) are then expanded in the 
normal modes of the system, the expansion coefficients 
become the photon annihilation and creation operators. 

In the Keldysh field-theoretical approach [33, 34] we 
calculate the partition function 



Z = TV([/p), 

where p is the density matrix at time t 
Trp = 1, and 



(5) 



-oo, with 



[/ = Tc exp 



-I I H{t)dt 



h=i. 



(6) 



is the time-evolution operator along the Keldysh con- 
tour C. The contour begins at t = — oo, where the state 



of the system is known, then goes forward in time up 
to < = oo, where it turns back and goes to t = — oo. 
Tc denotes the time ordering along the contour. The 
Hamiltonian H{t) is switched on adiabatically, starting 
from a trivial Hamiltonian H{—oo). H{t) is the same 
on the forward and backward branches of the contour. 
This condition leads to i7 = 1, and, hence, Z = 1. If 
the source terms that are different on the two branches 
are added to the Hamiltonian then Z ^ 1. The (func- 
tional) derivatives of the type 5Z[J]/SJ\j=o with respect 
to the sources J generate averages with the density ma- 
trix propagated from t = — oo to relevant times. 

The partition function can be written in the form of a 
functional integral over the fields (i.e., the classical func- 
tions) A and A^. To this end, we represent the classical 
action along the Keldysh contour as 

Sc = T^ [Al G-'A+- Al G-'A_] 
Ion 

= _i_[(^ci)t Gf-i^q + (^q)t G-^A% (7) 

where the subscripts "±" denote the fields on the forward 

and backward branches of the contour and the so called 
classical and quantum fields are defined by 



A'> = -^(M-A_). (8) 



The minus sign in front of the A[_ G~^A- term in Eq. (7) 
takes care of the time reversal on the backward branch, 
whereas A_ is the representation-free (vector) notation 
for the function A_{r,t) with the forward time ordering. 
It is convenient to consider A'^^ and as components of 
a single field 



A = 



(9) 



in the Keldysh space, which is twice the size of the origi- 
nal Hilbert space. (We will furnish the vectors and oper- 
ators in this space with a hat.) Then the contour action 
can be written in the form (dropping the subscript "C" ) 



iOTT 



(10) 



where G~^ has a 2 x 2 matrix structure with zeros on the 
diagonal and equal off-diagonal elements. 

In order to use S[A, A^ in the functional integral for Z, 
the operator G~^ has to be regularized [33] by imposing 
the causality structure on its matrix: 



G 



(G-i)A 
(G-i)i^ (G-i)K 

(G-i)^'A(r) = e(r, to) co^ + V^± iO+, 

(G-i)K = (G-^)^F - F(G-i)^. 



(11) 

(12) 
(13) 



Here, (G~^)^'^'^ are the retarded, advanced, and 

Keldysh components of the inverse Green function opera- 
tor. The operator F that parameterizes (G~^)^ depends 
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on the thermal distribution. Equation (12) is written 
under assumption of real e(r, w). In the medium with 
absorption, it is generalized to 

{G-')^^\r) = e'{r,uj)oj^ + ± ie"{r,co)uj^ (14) 
e'(r,a;)=Re[e(r,a;)], e"(r, w) = Im [e(r, w)], (15) 

where e" > 0. In the case of the gain medium, e" < 0, the 

time-dcpcndcnt Green function is exponentially diverg- 
ing, and its Fourier transform to the frequency domain 
does not exist. The frequency representation can be de- 
fined with the help of the Laplace transform, which is 
equivalent to introducing a fictitious absorption to the 
system. When the results of a calculation are trans- 
formed back to the time representation, they should not 
depend on the fictitious absorption. This means that one 
can perform the calculations in the frequency domain as- 
suming e" > 0, and obtain the final results by analytic 
continuation to e" < 0. 

The functional-integral representation of the partition 
function becomes 



Z =J^ I D\A,A^e'^^^'^'\ (16) 



where J\f is the nonessential normalization constant that 
ensures Z = 1 and the measure is defined by 

D[A,A^] = n 4Re^^.(r)) d{lmAi{r)) 

r.iu. 



Here and below we sot to unity the step size for the grid 
used to discrctizc the integral. 

We note that the Keldysh formalism is especially ap- 
propriate for the description of systems with absorp- 
tion or gain as it naturally takes into account the finite 
e"(r,w) in the causality structure of the inverse Green 
function, Eqs. (11), (14), and (13). 



B. Disorder average 

We shall study the effect of disorder in the refractive 
index and assume the absorption or gain in the system to 
be spatially rmiform. Hence, we represent the dielectric 
constant in the form 



e(r,w) = e'(w) -|- Ae'(r,w) + ie"{uj) 



(18) 



with the averages over disorder realizations 
(Ae'(r,w)) = and (Ae'(r, w) Ae'(r', w)) cx 5{v-y'). 

In order to define the scattering time r, let us, for 
a moment, neglect e"(w). The wave equation (1) can 
be interpreted as a time-independent Schrodinger equa- 
tion with the energy E{ijj) = t'{ixi)uP' and the potential 
energy F(r,w) = — Ae'(r, w) w^. The scattering time t 
and other characteristic time scales of the system are 
assumed to be much larger than liJ^^ ^ where cjo is the 
typical optical frequency. The slowly varying amplitude 



A(v, t) = A{r, t) exp{iLJot) satisfy the approximate equa- 
tion 

~8A ~ 
in— = [-V^ + V{v, ujo) - £{uJo)]A, (19) 

which is the time-dependent Schrodinger equation with 
the "optical Planck constant" 



~ _ d£{ujo) 
du){) 



(20) 



The scattering time can now be defined by analogy with 
the quantum scattering time via the correlation func- 
tion [35] 



(y(r,a;o)V(r',a;o)) 



2'KVT 



5{v-v'), (21) 



where v = dn/d£ is the quantum density of states (per 
unit volume). Note that u = uq/H, where uq = dn/duj is 
the standard optical density of modes. 

A disorder average of the partition function can be 
obtained by evaluating the functional integral 



{Z} = J D[V]Zexp J drV^{r,cuo) 

D[V] = l[^dV{r,uJo). 

(17) 

The disorder-dependent part of the action is 

AS[A,Aiv]=-±-A^v^A, 7-(;;) 



(22) 
(23) 

(24) 



where V{r,ujo) appears as an operator V diagonal in r; 
F is assumed to bo diagonal in r, as well. By completing 
the square, we obtain the disorder contribution to the 
partition function 



(e*^^) = exp 



AlTVT 



. (25) 



The short-hand notation A{y) is used for the Keldysh- 
space vector A with the fixed index r, i.e., it is a vector in 
the space with the reduced dimensionality; in this nota- 
tion, (v)A{y) involves a summation over the remaining 
indices, e.g., u and the Keldysh index. 

The negative sign in the exponent (25) is essential for 
the properties of nonlinear sigma model in the optical 
medium. In contrast to a fermionic system, the sign can- 
not be changed by commuting the fields. 



III. NONLINEAR SIGMA MODEL 

A. Hubbard-Stratonovich transformation. Saddle 
point 

The term of the fourth-order in the fields in Eq. (25) 
can be converted to a second-order term with the help of 



4 



the Hubbard-Stratonovich transformation yielding 



exp 



(l67r' 



j D[Q]e 



(26) 
(27) 



i=cl,q' 



The auxihary field Q is the Hermitian operator diagonal 
in r. The measure D[Q] is defined over the independent 
matrix elements by analogy to Eq. (17). The normaliza- 
tion constant Mq is determined by setting A = and 
A' = 0. The negative coefficient in front of TrQ^ de- 
termines the scale of Q and can be chosen freely. The 
present choice leads to the simple form of matrix A in- 
troduced in Eq. (34). In order to prove Eq. (26) more 
easily, one can define the matrix .4(r) = A{r) ® A^r) 
where the tensor product applies to the Keldysh and co 
subspaces. Then one represents 

A^QA = Tr{QA^). (28) 

Now the field Q can be integrated out after completing 
the square. 

Using Eqs. (25) and (26) in Eq. (16) and integrating 
out the fields A and A^, we obtain the disorder-averaged 
partition function 



iS[Q] 



-Tr 



(29) 
(30) 



where Gq is the inverse Green function operator that 
does not include the; disordered part of the dielectric con- 
stant and all Q-independent factors are included in the 
normalization constant A/g. 

In the limit of large scattering time, the main contri- 
bution to {Z) comes from the neighborhood of a saddle 
point. The saddle-point equation 



(31) 



follows from the condition of stationary variation of S[Q] 
with respect to ^Q. In the (k, co) representation, 

(Go ')S'^(k) = £{co) -k^± ie"{uj) (32) 

is diagonal. We will look for the solutions Q^'^ in the 

(cl, cl) and (q, q) blocks of Q, respectively, which are uni- 
form in r and diagonal in w. Equation (31) yields for 
these blocks: 



R.A 



(33) 



The sum over the modes can be converted into an in- 
tegral over i/dE, where £ = k^. In the limit uit ^ 1 
and e" <^ e', the lower integration limit can be extended 
to — oo. Then {Q^''^)uj = ±« is the solution. The full 
matrix can be written in the form 



Qo = iA, A 



1^^ 2F 
-1^ 



(34) 



which includes the regularization in 1^'^ = e='=*"'^^ and 
the Keldysh block. The regularization leads to an impor- 
tant property Tr = 0. 

We note that the saddle point Qo lies outside of the 
manifold of Hermitian matrices Q. The diagonal part 
of Qo is anti- Hermitian; this property can be traced back 
to the negative sign in the exponent in Eq. (25). The 
Q manifold can be continuously deformed to pass through 
the saddle point by making the transformation Q i— > e^'^Q 
in the neighborhood oi Q = A. An (f) changes from 
to 7r/2, the logarithm argument in Eq. (30) has no zero 
eigenvalues if e" > 0. Thus, no singularities are crossed 
by ex.p{iS[Q]) during the deformation. 



B. Effective action 

The main contribution to (Z) arises from the fluctu- 
ations about the saddle point that satisfy the condition 
TrQ'^ = 0. Such fluctuations produce weak variations of 
the action S[Q] (30). The matrices Q having the above 
property can be represented in the general form 



Q{r)=iR-\r)AR{r), 



(35) 



where R is diagonal in the r representation. 

In what follows we present the results of the calculation 
and refer the reader to the Appendix for details. After 

substituting the parameterization (35) in Eq. (30) and 
omitting the Q-independent contribution, we arrive at 

iS[Q] = -Tr In (i + g^R[^Go\ R-^]^ 

^ -Tr {g^R[^Go\R-']) 

^Tr {gjR[jG^\ R-']y . (36) 



+ 



where 



g=[G^' + t^iA 



(37) 



is the disorder-dependent Green function operator [see 
Sec. IV B]. The action is expanded in the fluctuations 
about the saddle point, which are described by the com- 
mutator [^Gq^ , R~^]; at the saddle point ^ = i the 
commutator vanishes. The disorder-free inverse Green 
function consists of the conservative and nonconservative 
parts: 



7(Goi)^(k) = [£{oj)-k^]i + ie"{u)oj^A. 



(38) 
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The latter results in a nontrivial contribution to the com- 
mutator due to the Keldysh structure of A. 

There are three leading-order contributions to S[Q]. 
Using the £{u}) part of ^Gq^ in the linear term in Eq. (36) 

we arrive at 

iSi[Q]^tTTiyhTr{dtQ), {dtQ)u ^ {dtQtt')t'=f (39) 

The contribution of the k"^ part of jG^^ to the linear term 
of Eq. (36) is neglected compared to its contribution to 
the second-order term, which gives 



iS2[Q]c^--vhDTr{drQy. (40) 
To derive this result, we used the property [35] 

h 



2ttv 



EC(k)C(k) 



1 2e"wg 



(41) 



This condition specifics the regime when the light propa- 
gation is diffusive. If this requirement is not fulfilled, the 
massive fluctuations beyond the NLSM have to be taken 
into account. 



IV. LIGHT DIFFUSION 

In this section we calculate the disorder-averaged 
Green-function correlator. In particular, we consider 
the contribution that arises from the fluctuations of the 
field Q in the neighborhood of the saddle point. The cor- 
relator possesses a diffusion-pole structure modified by a 
finite absorption/gain rate. 



A. Fluctuations about the saddle point 



yielding the effective scattering time r and defined the 
effective diffusion coefficient in two dimensions. 



D 



1 



,2 



(42) 



where v is the group velocity of light in the medium. In 
Sec. IV we show that f and D are the relevant parameters 
to describe the diffusive light propagation [see Eq. (75)]. 
Finally, the nonconservative part of 7G0 ^ yields, in the 
linear order in the commutator, 



iSslQ] = ■Kve"ujl Tr {ikQ + A^). 



(43) 



The contributions 51^2,3 [Q] sum up to yield the NLSM 
effective action 



iS[Q] = 
- TTZ^oTr 



II, ,2 



e w. 



■(iAQ + A^) 



(44) 



The action vanishes at the saddle point, S[iK\ = 0. The 
key assumption behind the NLSM is the smallness of the 
action for fluctuations of Q restricted to the manifold 
TrQ'^ = 0, compared to the action for arbitrary fluctua- 
tions about the saddle point. The terms S'i_2[Q], which 
also appear in the NLSM for disordered fermionic sys- 
tems [33], depend only on the derivatives of Q. There- 
fore, the dominant contribution to the partition func- 
tion comes from the fluctuations Qtt' (r) [or Rtf (r)] that 
are slowly varying functions of r and {t + t')/2. These 
"massless modes" are associated with the diffusive light 
propagation. The assumption of slow variation justifies 
neglecting of higher-order terms in the expansion (36). 
The contribution 6*3 [Q] results from the nonconservative 
nature of the medium. It is, in general, comparable to 
the to the "massive" Tr term, unless the rate of ab- 
sorption or gain is smaller than the scattering rate: 



(45) 



We consider the parameterization 



Q = i/7e-^/2a,e^/2(7-i, U = (j- 



1 F 
-1 

where ct^ is the Pauli matrix. Because A = UdzU~^ (if 
the regularization of unit operators is neglected), this 

parameterization is equivalent to Eq. (35) with R = 
?7exp(VF/2) J7^^. It can be verified by explicit calcu- 
lation that the diagonal blocks of W do not contribute 
to S[Q] and the Green- function correlator, at least, up 
to the second order in W. We, therefore, represent this 
field in the form 



W 



w 
t 



w 



(47) 



The specific choice of as a Hermitian matrix is justified 
by the requirement of convergence of the functional inte- 
gral for the partition function (see below). The operator 
w is diagonal in the r representation. 

By expanding the parameterization (46) in the powers 
of W we find the first- and second-order deviations from 
the saddle point, 



5Q 



(1) 



-w — Fw^F 
w'^F 

i ( wvj^ ww'^ F + Fww^ 



w 







r,t 



(48) 
(49) 



We note that only the latter matrix has the causality 
structure; however, the fluctuations of Q are not required 

to obey causality. By using (JQ^^'^^ in Eq. (44) we can 
calculate fluctuations of the effective action. 

The flrst-order variation of S[Q] depends on the deriva- 
tives [36] of the distribution function F generated by the 
flrst two terms in Eq. (44) ; the third term yields an identi- 
cally vanishing flrst-order contribution. The saddle-point 
equation (33) determines the retarded and advanced sec- 
tors of Q, but not the function F. By setting to zero the 
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variation of the effective action near the saddle point, we 
obtain the Usadel equation 



(50) 



for Fw (r) in the mixed representation of the slow time 
variable t = {t + 1')/2 and the large frequency u, conju- 
gate to t ~ t' . 

The second-order variation is 



-TTZ/oTr 



D 



2 h 



2 



E 



'(k)l' 



-z(u;-a;')+5fc2 + ^i^ 



(51) 



where w;(k) is the Fourier transform of w{r). Of the 
two terms with spatial gradients, the second term has a 
zero trace. The first term yields the Dk^ contribution to 
SS^'^^ , as well as the additional correction 



.^Dl,[w\d.F)]\ (52) 



where "tr" denotes the trace of operators that do not 
have the Keldysh matrix structure. This correction van- 
ishes when F{r) is uniform, which we will assume. With 
the help of Eq. (51), the disorder-averaged partition func- 
tion can be approximated by the functional integral 



{Z)^N^ j D[w,w^]e 



tl JSS^^\w,w^] 



(53) 



where Nw is a normalization constant. For a medium 
with gain, the divergence of the integral for the modes 
with 



k<k 

n 



2e"ujl 



hD 



(54) 



indicates that the long-scale fluctuations become unsta- 
ble due to onset of lasing (see Sec. IV C). Thus, in 
the long-wavelength limit the linear-gain theory breaks 
down and the nonlinear effects have to be taken into ac- 
count [15, 16]. 



with respect to the source fields: 



(52Z[J, Jt] 



Qjk ^ 2) g'™ (3, 4) G^" (1,4) G"= (3, 2) 
1 6^Z[J, Jt] 



(55) 



(167r)2 S[Ji{l)]* (5Jfc(2) 6[Ji{3)]* (5J™(4) 
Z[J, P\=N j D[A, it] e'S[^.^']+-^'^+^'^, (57) 



J=jt=o 
(56) 



where j,k,l,m = cl, q and 1,2,... are full sets of co- 
ordinates in some representation, e.g., 1 = (ki,wi), 
etc. By inverting the matrix (11), we identify the sec- 
tors of the Green function as C^''^ = G^, C^'"^ = G^, 

The disorder-averaged Green functions and correlators 
are obtained by using the above expressions with the 
disorder-averaged partition function [37] 



{Z[J,P]}=Mq Id[Q] exp(i5[0] + 167rJt^^jj, 



We find, in particular 

iR,A,K 



(g^'^'^(i,2)) = (c;J^'^(i,2))q, 

(G^(l,2)GA(3,4)) 

= (5^(1,2)5^(3, 4)-fe|(l,4)e^(3,2))Q, 



(58) 

(59) 

(60) 
(61) 



where the average {■■■)q over Q is performed with the 
exponential weight exp(zS'[Q]). Equation (60) shows that 
Q — Q^^ [Eq. (37)] is the disorder-averaged Green func- 
tion in the lowest-order saddle-point approximation. The 
component = Gq"^ in Eq. (61) is, in general, non-zero 

when Q does not have the causality structure. This ob- 
servation is essential for the following calculation. 
We calculate the Green-function correlator 

7e(l,2,3,4) = (GI^(1,2)GA(3,4)) - {G''{1,2)){G^{3,4)) 

(62) 

by expansion about the saddle point. The lowest-order 
correction to the Green function (59) is 



B. Disorder-averaged correlator 

Green functions and their combinations can be ex- 
pressed in terms of derivatives of the partition function 



2r 



Gj6Q^^'>G 



2t 



G^wG^ + FG^w'fG^F FG^w^^G^ 
-G^w^G^F -G^w'^G^ 



■ (63) 
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The Gaussian averages with the action (51) are as follows: 

{w)^={w%=0, («;t(l,2)«;t(3,4))^ = (64) 
(w;(l,2)w;t(3,4))^ 

2 ^ki — k2 ,k4 — k3 ^UJi ,UJ4 ^UJ2 ,^3 



^'^0 D(ki - k2)2 - ^(wi - C02) + ^ ' 



(65) 



Therefore, the leading contribution to the correlator 
comes from the K-Q term in Eq. (61), which is given 
by the product of diagonal blocks in Eq. (63). We find 



n{i, 2, 3, 4) = ^—-^g^{i) g^{2) 0^(3) g^{A) 

^ '^ki— k4,k2-k3'^a;i,ui2^t>'4,a'3 

5(ki - k4)2 - i{cOi - COi) + ^ ■ 



(66) 



The correlator has a diffusion pole with the diffusion co- 
efficient D. The pole is modified by the e" term that 
arises from the corresponding contribution in the effec- 
tive action (44). This term defines the absorption rate 



1 _ 2e"a;2 



(67) 



negative for gain. 



Discussion 



The pole structure of the correlator implies that the 
light intensity / in the medium satisfies the diffusion 
equation with a nonconservative term: 



{dt-DV^ + T-')l = 0, 

- 1 . 1 9 /I 1 

2 2 \t 

We compare this equation with 



D' = -v^ 
2 



1 2 
- + — 



(68) 
(69) 

1 = 0, 
(70) 

(71) 



that follows from the photon transport equation (see 
Eq. (15) of Ref. [17]). According to Ref. [17], the light 
propagation is diffusive if the second derivative with re- 
spect to time in Eq. (70) can be neglected. This is the 
case when 



T < At, 



(72) 



where At is the characteristic time scale of intensity vari- 
ation. The reaction of the medium on a fluctuation of 
intensity will be determined by the shortest time scale. 



so that At < |Ta| can be assumed. Therefore, when ne- 
glecting the corrections of the order of r/At in Eq. (70), 
we also have to neglect the contributions of the order 
of r/ra. In particular, it is consistent with the diffusion 
approximation to set 



(73) 



The independence of absorption for the diffusion coeffi- 
cient was also supported by the numerical evidence in 
Ref. [17]. It is worth commenting on the claim [18, 19] 
that the diffusion coefficient in the medium with absorp- 
tion must be equal to D even for r/ra ~ 1. A closer 
look at the derivation of the diffusion coefficient from 
the transport equation in Ref. [18] reveals that the time- 
derivative terms neglected in Eqs. (A9) and (All) of that 
article would yield the diffusion coefficient 



D" = -V^T 

2 



1-2- 



D' 



(74) 



were they taken into account. Thus, the (approximate) 

independence of the diffusion coefficient of absorption is 
a consequence of the self-consistent application of the 
diffusion-approximation conditions (45) and (72). 

The NLSM effective action (44) is derived under the 
condition (72) as well. This condition guarantees the slow 
variation of Q, and makes it possible to neglect the con- 
tribution of £{uj) part of jGq^ [Eq. (38)] to the second- 
order term in Eq. (36). This contribution would result 
in a second-time-derivative term in the effective action. 
Again, the diffusion approximation requires that we set 



(75) 



in the NLSM expressions. Thus, the NLSM and the the- 
ory of transport equation agree in the diffusive regime. 

In the medium with gain, the diffusive relaxation com- 
petes with the amplification. Because the long-scale in- 
tensity fluctuations disperse slower, they become unsta- 
ble, and the random lasing sets in. The cutoff wavenum- 
ber kmin (54) determines the critical sample size 



i = Vd'\ 



(76) 



above which the system is lasing and the linear-gain the- 
ory does not apply. Alternatively, the above expression 
yields the lasing-threshold value of [raj if / is given. 



V. CONCLUSIONS 

We obtained the functional-integral form of the par- 
tition function for an optical medium with linear ab- 
sorption or gain. Keldysh technique is particularly suit- 
able for description of nonconservative systems because 
it provides a natural representation for the action. The 
disorder- averaged partition function is expressed as a 
functional integral over the auxiliary matrix field Q. 
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Within the framework of nonhnear sigma model, we con- 
sidered the fluctuations about the saddle point that fulfill 
the condition Tr — 0. We found that the effective ac- 
tion ^[(5] for these fluctuations contains an extra term 
due to absorption or gain. 

With the help of the nonlinear-sigma-modcl parti- 
tion function, we computed the disorder-averaged Green- 
function correlator. The leading contribution from the 
vicinity of the saddle point has the diffusion-pole struc- 
ture modified by a finite absorption/gain rate. The diffu- 
sion coefficient is found to be approximately independent 
of the absorption or gain in agreement with the theory 
of photon transport equation. In the medium with gain, 
the linear theory is not applicable in the long- wavelength 
limit. If the sample size exceeds a certain critical length, 
the random lasing sets in. 



follows from the saddle-point condition; furthermore, 
X^k (k) k = due to the symmetry. After calculat- 
ing the Ak sum we arrive at 



q2l 6-1 



= TT!/ Tr j dr dt dt' iARftir) [ihdt + d^] R'} (r) . (A.4) 

Applying the representation (35) we obtain Eq. (39) from 

the dt part. The 9^ part is neglected compared to iS'2[Q]; 
the latter contribution is multiplied by HD ~ coof ^ 1. 



Appendix: Derivation of the effective action 
1. Derivation of Eq. (36) 



3. Derivation of Eq. (40) 

We substitute the part of 7G0 ^ in the second trace 
in Eq. (36) to get 



After substitution of the parameterization (35) in 
Eq. (30) we obtain 



iS[Q] = -Tr In 



R^R-^ ( R^Gq^R-^ + i—A 



= -Trln7-Trln(7^-i+i?[7Go\i?-i]) . 

(A.l) 

By separating Tr ln(7^r~^) and dropping the Q- 
independent terms we arrive at Eq. (36). 



2. Derivation of Eq. (39) 

The conservative part of jGq^ , when substituted in 
the first trace in Eq. (36), yields 



iS2[Q] = ^Tr ^ ^(ki)7^(ki - \^2)R-\^2 - ks) 



ki...k4 



X a(k3)7i?(k3 - k4) i?-l(k4 - k,)ik'^ - ki){ki - fcf) 

~2Tr ^ g(i)^RfAh-^^^^\R-\Ak') 

1~ A 1~ ^ 



kAk 
AkAk' 



X miR{ -Ak - ^]i±^ ) ^-i(Ak) 



X (k- Ak)(k- Ak'), 



(A.5) 



where k = ki/4 and we neglected the contributions 

of higher order in Ak = (ki + k4 — k2 — k3)/2, Ak = 
k4 — ki, and Ak' = k2 — ks. We use the representation 



g^ = ^g^ii + A) + i{i-A)g^ 



(A.6) 



iS^lQ] =-TrJ2 ^c.(k)7^c.c.'(k - k') ^-/„(k' - k) 



k,k' 



X [£{uj') - £{w) - k'^ + 
-Tr J2 4(k)7^c.,c.+Ac.(-Ak) RzIao^J^^) 



k,Ak 



X [/lAw- (2k-F Ak) • Ak], 



and the well-known relations [see Eqs. (42) and (41)] 

E - l^^'f^DS^e^ (A-7) 
k 

E a^(^) (k) (k) k^k^ ~ 0, (A.8) 



(A.2) to find 



where Aw = 



and Ak = k' - k. We note that R 



is peaked at small wave vectors in the k representation. 
The sum 



(A.3) 



iS2[Q] = --TTiyhDTT[{l + A)R{drR'^) 
■ (i - A)^ (drR-^)] = ^TTvllD Tr [dr{R-^AR)]^, (A.9) 

from which Eq. (40) follows. 
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4. Derivation of Eq. (43) 

The nonconservative part of jGq^, being substituted 
in the first trace in Eq. (36) yields 

iSalQ] ^ - ie"u>l Tr ^ ^(k)7i?,(k - k') 

kk' 

X [A^-^(k' - k) - R-^{k' - k) A]. (A. 10) 

Wc change the variable k' = k + Ak and apply Eq. (A. 3). 
After cyclically moving the operators under the trace we 
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